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Abstract 

We evaluate lower bounds on the sum of the up and down quark masses. The bounds follow 
from the constraints provided by the dispersion relation obeyed by the two-point function of the 
scalar current density when combined with properties of the scalar spectral function at long- 
distances and perturbative QCD at short-distances. Our results point to values of niu + rrid 
somewhat higher than those reported recently in the literature using lattice QCD simulations. 
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I Introduction 



The light u,d and s quark masses are fundamental parameters of the Standard Model of Particle 
Physics whose precise values are still affected by large uncertainties. The difficulty arises from the 
fact that it requires to have a good handle on QCD non-perturbative effects in order to extract their 
values (much smaller than the typical hadronic scale of a few hundred MeV) from the properties of 
the hadrons in which they are confined. The masses we refer to are the ones in the QCD Lagrangian 
and it is customary to report their values in a mass-independent subtraction scheme such as the 
MS-scheme at a reference scale /i which has conventionally been fixed by the lattice community at 
/i = 2 GeV, and which we also adopt. 

The extraction of the light quark masses from physical observables has been attempted using 
three independent methods: spectral function sum rules for hadronic correlation functions, lattice 
simulations, and chiral perturbation theory (ChPT). A summary of the spread of results obtained 
(also within each method) can be found in the successive versions of the Particle Data Book [Ij with 
references therein. The suggestion that ChPT can be used to extract ratios of the light quark masses 
from the physical properties of the low pseudoscalar particles goes back to earlier work by Weinberg [2] 
which has been subsequently improved (see e.g. ref. [3] and references therein). ChPT can thus relate 
determinations of m„ + md to those of m„ + nis in a reliable way. 

We shall be concerned with the question: how small can the light quark masses be? This question 
was already addressed in an earlier paper ^4. where a variety of lower bounds were derived. Here 
we propose to reexamine the bounds concerning the scalar-isoscalar channel where a number of new 
developments have appeared in the meantime. On the one hand, from a phenomenological point of 
view, the scalar pion form factor which governs the dominant two-pion contribution to the spectral 
function in the scalar two-point function, has been the object of careful analyses (see section III 
below). On the other hand, from a theoretical point of view, the two-point function in question has 
been evaluated at the five loop level f5] ^51 in perturbative QCD (pQCD), thus providing a better 
handle on the choice of the euclidean Q^-value at which one can be confident that the perturbative 
regime applies. 



II Bounds from the Scalar Channel 

The scalar density current operator in question is 

S{x) — rh [uu + dd] (x), rh= — 
with 

*(g2) = 1 fd^x e'^-^0\T{S{x)S{0))\0) 



(2.1) 



(2.2) 

the associated two-point function. This function obeys a dispersion relation which in QCD requires 
two subtractions. In terms of its second derivative, which gets rid of the two subtractions, and for 
Euclidean values of = —q^, the dispersion relation then reads 



:Re^{q^) 



dt— 



-Im*(t) . 



(2.3) 



This is the master equation from which the bounds on rh will be extracted. 

The behaviour of the '^"{Q^) function in pQCD at large euclidean values: » Aq^^, is known 
from more and more refined calculations (see refs. [3 to [ll] and [SI [6]) 



*"(Q') = 



Nc rn^{Q^) 
47r2 Q2 



2^ r 



llas(g') 



(2.4) 



where the dots represent higher order terms which are now known to 0{a\). Similarly, the spectral 
function ilm^'(t) on the r.h.s. of Eq. (|2.3p is also known in pQCD by analytical continuation of ^'(g^) 



1 



to large time- like values of = < >> ^qcd order ©(af) [S], 



^Ini*(i)|pQCD = ^ rn^it) t 



In terms of physical degrees of freedom, the spectral function is given by the sum 



(2.5) 



(2.6) 



where the sum over T extends to all possible on-shell states of the hadronic spectrum with the quantum 
numbers of the scalar current, including the integration over their phase space. These contributions 
add positively and, therefore, if the sum is restricted to the lowest possible contribution of two pions 
only, there will follow a rigorous lower bound to the full spectral function and hence, via the dispersion 
integral in Eq. (j2.3p and the definition of the S{x) current in Eq. (|2.ip . to the quark masses. More 
precisely 

11 3 
Im^'(t) > -Im*(t)|r=7r 
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where F{t) denotes the J = 0, / = scalar pion form factor defined by 

{iT%p)n\p')\s{m = ^'''m, t = {p+p'f. 

Therefore, for sufficiently large values of where pQCD in the evaluation of '^"{Q^) applies, 
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As discussed in ref. [T, knowledge of the scalar pion form factor at the origin F{Q) and of its mean 



squared radius 



m = F{(i) 



(2.10) 



results in a lower bound for as a function of Q^. A more restrictive bound was also derived in [4] 
using the fact that the phase of the scalar form factor Spit) in the elastic region 4m^ < i < 16m^ 
is precisely the J = 0, / = 07r — tt phase-shift 5^{t). Variations on similar analyticity properties 
concerning the strangeness changing scalar form factor have since then been also discussed in the 
literature (see e.g. ref. [12] and references therein). 

The replacement of the full scalar spectral function by the r.h.s. of the inequality in Eq. (j2.7p 
is expected to be a gross underestimate at large values of t. Indeed, for t large, pQCD predicts 
the spectral function to grow as 0{rh?{t)t), while the scalar pion form factor drops as 0{\/t) up to 
logarithms [131 114] . This suggests an improvement of the inequality in Eq. (|2.9p in the following way. 
First we separate the dispersion integral in the r.h.s. of Eq. (|2.3p in two pieces 



J(2m„)2 



(i + Q2)3 ^ 



-Im*(t) 
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[t + Q2)3 



(2.11) 



with the scale io chosen sufficiently large so that the spectral function in the high energy integral 
to t "^^ oo can well be approximated by its pQCD expression. Then we use the spectral function 
inequality (j2.7p only in the low energy integral (2to^)^ < ^ < ^o- This results in the following improved 
inequality 
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We wish to emphasize that this inequahty is rigorous, provided that both and to are chosen suf- 
ficiently large, so that pQCD can be applied to the evaluation of the function in the deep 
euclidean as well as to the dispersion integral in the tg < t < oo interval. We next proceed to the 
evaluation of this bound. 

Ill Evaluation of the Improved Bound 

The basic ingredient in the evaluation of the low-energy integral in Eq. (|2.12p is the Mushkhclishvili- 
Omnes representation of the form factor 

|F(i)| = |f(0)|expj-P r ds^^] , (3.1) 

[ 74„2 S(S - t) J 

which follows from the analyticity properties of F{t) if the form factor has no zeros and goes as 0(l/t) 
(up to logarithms) at large t. Then the phase (5f(0 must approach tt asymptotically!^. The form 
factor at the origin is known from ChPT [16] 

F(0) = ml 

where the error here takes into account a generous error on the low-energy constant Z3 as well as on 
higher order corrections [17]. As already mentioned, the phase Sp{t) in the region 4m^ < t < 16toJ is 

3 . 5r 

3 ; 
2.5; 

T3 9 - 
U 

[n 1 • 5 " 

1 ; 
0.5; 

0.5 1 1.5 2 

[Gev] 



327r2F2 



(0.99 ±0.02), 



(3.2) 




Figure 1: Shape of the phase SF{t) as a of function oft. In the interpolation region the error varies 
smoothly up to 20% at Vt ^ 1.7 MeV. The errors below the K — K threshold are of the same size as 
those given in ref. \15^ . 

given by the tttt phase-shift 5o(i) and it has been extracted from tt— tt scattering experiments |18j . from 
2 For a detailed discussion see e.g. ref. I15| . 
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Ki4 decays [12] and K ^ tttt decays 20 . In fact the relation Spit) = i5o(i) can be extended up to the 
KK threshold, since inelastic processes do not play a major role in the interval 16m^ < t < 4M|. |21j . 
Excellent fits to the phase Spit) in this region can be found e.g. in refs. The opening of the KK 

threshold, however, produces a square root singularity at t — 4M|- which causes a dip in the elasticity 
in the region 1 GeV < Vt < 1.1 GeV. The corresponding effect in the determination of the phase Spit) 
which appears in Eq. (j3.ip was first discussed in ref. [24] within a two-channel representation, where 
only the tttt KK transition is assumed to drive the whole effect of the inelasticity. It was later 
confirmed by a more elaborate analysis in ref. |25j , where a variety of parameterizations [26j are taken 
into account. As a result, the overall shape of the phase Spit) is now rather well known up to energies 
Vi < 1.4 GeV. Beyond that, it can be assumed to join smoothly the asymptotic regime predicted by 
pQCD [131 E]. In Fig. 1 we show the plot of the phase Sp{t) thus obtained (the solid line), where 
we also show a cautious error margin of ~ 20% in the interpolation region: 1.4 GeV ^ ^ 2.2 GeV 
(the dotted lines) . Fortunately, this alleged source of uncertainty turns out to have little influence on 
the light quark mass bounds, as far as the interpolation is smooth. This is because the low-energy 
integrand in Eq. p.l2p is modulated by three powers of {t + Q'^) in the denominator which, because 
of the large Q^-values we are considering, makes the integral rather insensitive to the precise values 
of the phase in this intermediate region 
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Figure 2: Shape of the low-energy integrand piow {Q^,t^ (the solid curve) in Eq. i3.4\ ) and the high- 
energy integrand phigh {Q^,t) (the dotted curves) in Eq. i3.5\) as a function of the t for the choice 
^4: GeV^ . The four dotted curves correspond to the successive pQCD approximations. 



^We wish to recall that the value of Sp(t) in this intermediate region has been a source of debate in the literature 1151 
I27II28II29] . where it is argued that beyond \/i — 1-4 GeV other channels like the four pion channel may sensibly distort 
the shape of the phase. The debate focuses on the precision that the mean squared radius (r^)J can be extracted from 
the integral 

\F{0)\ -IT J Ami ^ 

We would like to stress the fact that, contrary to this observable, the integral we are concerned with is much less 
sensitive to extrapolation uncertainties. In any case, the result we obtain for the mean square radius, using the solid 
curve in our Fig. 1 is in good agreement with the value = (0.61 it 0.04) fm^ quoted in |21| . 
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The solid curve in Fig. 2 shows the shape (without errors) of the low-energy integrand in Eq. (|2.12p 

Plow {Q'^,t) 



(3.4) 



as a function of t for a reference choice = 4 GeV^, with the form factor \F{t) \ evaluated as described 
in the previous paragraph. The dashed curves in the same figure show the shape of the high-energy 
integrand: 



Phigh {Q^,t) = 



4 m^jt) t 
3 ^(t + Q2)3 



(35.94- 1.359 n,)^"^(*) 



as{t) 



(164.14 - 25.77 nj + 0.259 n) 
(39.34 - 220.9 n/ + 9.685 n) - 0.0205 n'j) 



as(t) 



O 



as{t) 



(3.5) 



at the same reference choice = 4 GeV^, for different approximations of the pQCD series in powers 
of "'J*'' ; i.e., from one power (the curve with the largest dashing) to four powers (the curve with 
the shortest dashing). Each curve here is modulated by the appropriate running quark mass 'm?{t) 
resulting from the inequality in Eq. ()2.12p and approximated at the corresponding number of loops. 
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Figure 3: Plot of the lower hound (the gray hand) ohtained for [m,i + md\{^J. — 2 GeV) in MeV as a 
function of in GeV^, where we have also set to = in Eq. h2.12]) . The width of the band reflects 
the propagation of errors, mostly from the determination of the phase 6p(t) shown in Fig. 1. The 
horizontal dashed band corresponds to the most recent lattice results quoted in the text. 



We find that for t > 4 GeV^ the error between the ©(af) and the ©(a^) approximations to 
Phigh is less than 5%. It seems then natural to adopt the simplest choice to = for the 

separation scale in Eq. (|2.12p and proceed to the numerical evaluation of the bounds as a function 
of only one large scale. The corresponding lower bounds for [to„ + md\{^ = 2 GeV) as a function 
of Q2 ^ in the interval: 4 GeV^ < < 9 GeV^ are then shown in Fig. 3. The thickness of the 
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band indicates the effect of a scan in the propagation of the errors discussed above. They are largely 
dominated by the error in the extrapolation of the form factor. The choice > 4 GeV^ is already 
large enough for the error to be almost independent of Q^. In other words, with the input discussed 
above, the bound at = 4 GeV^ should already be a rigorous bound. 

IV Conclusions 

It is interesting to compare the values that we obtain for the lower bounds on m„ + irid, as shown in 
Fig. 3, to the latest determinations of the light quark masses reported from lattice simulations: 

mu + TOd = 6.6 ± 0.6 MeV MILC Collaboration [32] (4.1) 
mu+rrid ^ 6.4 ± 0.6 MeV HPQCD Collaboration [33] . (4.2) 

These numbers correspond to masses in the MS-renormalization scheme at a reference scale /j, = 
2 GeV, like our bounds in Fig. 3. The quoted error, is our addition in quadrature of the statistical, 
lattice systematics, pQCD, and electromagnetic/isospin effects quoted in the original papers. These 
two lattice results are in remarkable agreement with each other. They correspond to the dashed band 
in our Fig. 3. As seen from this comparison, they are somewhat smaller than the lower bounds. We 
find e.g., 

m„ + md > 9.7±0.4MeV at = 4 GeV^ (4.3) 

mu + md> 7.5 ± 0.3 MeV at = 6.5 GeV^ , (4.4) 

m„ + md > 6.1±0.3MeV at = 9 GeV^ (4.5) 

One has to go to values as high as g^ > 7 GeV^ to be in agreement. On the other hand, the bounds 
can very well accommodate some of the QCD sum rules results in the literature, like e.g., 

rriu + md = 9.2 ± 2.4 MeV [H] (4.6) 
m„ + md = 9.4± 1.8 MeV [35]. (4.7) 

Lower values than these ones, however, have also been reported more recently, like e.g., 

niu + rrid = 7.5 ± 0.7 MeV 36J . (4.8) 

It will be interesting to see in which direction the future determinations, both from lattice QCD and 
from QCD sum rules, will evolve. 
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